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Abstract 

We explore in more detail the modified quark meson coupling (MQMC) model 
in nuclear matter. Based on previous studies two different functional forms 
for the density dependence of the bag constant are discussed. For uniform 
matter distributions the MQMC model can be cast in a form identical to QHD 
by a redefinition of the sigma meson field. It is then clear that modifications 
similar to those introduced in QHD will permit the reproduction of all nuclear 
matter properties including the compressibility. After calibrating the model 
parameters at equilibrium nuclear matter density, the model and parameter 
dependence of the resulting equation of state is examined. Nucleon properties 
and scaling relations between the bag constant and the effective nucleon mass 
are discussed. 
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I. INTRODUCTION 


The description of the nuclear many-body problem in terms of strongly interacting quarks 
and gluons is one of the major challenges in nuclear physics. At present, however, rigor¬ 
ous studies of QCD are restricted to matter systems at high temperature and zero baryon 
density. Because of the nonperturbative features, it appears very difficult to derive from 
this theory predictions for processes at energy scales relevant for low- and medium-energy 
nuclear phenomenology. 

On the other hand, it is well known that despite these difficulties nuclear phenomenol¬ 
ogy can be efficiently described using hadronic degrees of freedom. While this framework 
has been very successful in describing the features of nuclear matter and the binding en¬ 
ergy systematics of finite nuclei, experiments, such as deep inelastic scattering off nucleons, 
provide evidence that the standard hadronic picture has to be corrected. For example, the 
prominent EMC effect which reveals medium modifications of the internal structure of the 
nucleon [|l|. 

Moreover, hadronic models are often extrapolated into regimes of high density and tem¬ 
perature to extract the nuclear equation of state, which is the basic ingredient in many 
astrophysical applications and in microscopic models of energetic nucleus-nucleus collisions. 
One can expect that under these extreme conditions quark degrees of freedom become im¬ 
portant. 

To address these issues it is necessary to build theories which incorporate quark-gluon 
degrees of freedom and which help to bridge the gap between nuclear phenomenology and 
the underlying physics of strong interactions. An important criteria for these new models is 
that they reproduce results based on the established hadronic framework. 

Guided by the symmetry breaking pattern of QCD, much effort has recently been devoted 
to the study of effective Lagrangians for low-energy strong interactions. A typical example 
is the Nambu and Jona-Lasinio (NJL) model [^. These models are mainly concerned with 
the spontaneous breaking and restoration of chiral symmetry but it is not clear if basic 
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features of nuclear phenomenology, such as saturation of nuclear matter, can be described 
properly. On the other hand, almost a decade ago, Guichon [Q proposed a quark-meson 
coupling (QMC) model in which nucleons arise as MIT bags interacting through meson mean 
helds. This model was rehned later by including center-of-mass corrections and applied 
to nuclear matter [p|-p!^ and also, more recently, to hnite nuclei pr3|JH 

Although it provides a simple and attractive framework to describe nuclear systems in 
terms of quark degrees of freedom, the QMC model has a serious shortcoming. It pre¬ 
dicts much smaller scalar and vector potentials than obtained in successful hadronic models 


11,12]. As a consequence the nucleon mass is much too high pl6, 11,12] and the spin-orbit 


force is too weak to explain spin-orbit splittings and spin observables in hnite nuclei. 

A well established framework for relativistic hadronic models is provided by quantum 
hadrodynamics (QHD) [|T^. Numerous calculations have established that relativistic mean- 
held models based on QHD provide a realistic description of the bulk properties of hnite 
nuclei and nuclear matter [|^]. One of the key observations in their success is that nucleon 
propagation in the nuclear medium is described by a Dirac equation featuring large scalar 
and vector potentials. They emerge from the interaction of a nucleon with all other nucleons 
in the Fermi sea via the exchange of isoscalar scalar and vector mesons. 

Recently, it was pointed out that the small vector and scalar potentials in the QMC 
model are due to the assumption that the bag constant does not change in the nuclear 
environment ]|TI| , p!2[] . By introducing a density dependent bag constant it was demonstrated 
that large scalar and vector potentials can be produced. A necessary condition is that the 
value of the bag constant in the nuclear environment signihcantly drops below its free-space 
value. As a consequence relativistic nuclear phenomenology can be recovered from a modihed 
quark-meson coupling (MQMC) model ]|TT| . 

The central issue of the MQMC model is the density dependence of the bag constant. A 
priori this is not known and the idea of the MQMC model is to parametrize the bag constant 
and to determine the parameters by calibrating to observed nuclear properties. Two different 
model types have been proposed |Tl|JI^: a direct coupling model in which the bag constant 
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is a function of the scalar field and a scaling model in which the bag constant is related to 
the effective nucleon mass. The density dependence is then generated self-consistently in 
terms of these in-medium quantities. However, the proposed models are not flexible enough 
to predict nuclear matter properties on a satisfactory level. The trend in the MQMC model 
is to predict reasonable values for the effective nucleon mass but compressibilities which are 
too high and vice versa 


We adopt the approach of Refs. [|^,0 with the goal of generalizing the proposed models. 
We assume two different functional forms for the bag constant. In one case it depends 
on the scalar held only; in our second model we study the bag constant as a function of 
the effective nucleon mass. To provide sufficient hexibility we model the functional form 
by using polynomial and Fade parametrizations. The unknown coefficients serve as our 
model parameters. We investigate how well the models can be calibrated by htting the 
parameters to properties of nuclear matter, which we take to be: the equilibrium density 
and binding energy (p^, —Cq), the nucleon effective (or Dirac) mass at equilibrium 
and the compression modulus {Kq). In general the models contain more parameters than 
there are normalization conditions. Thus families of models can be generated which describe 
exactly the same nuclear matter properties at equilibrium. This allows us to study different 
physical situations and to search for model dependence in the predictions. 

Because the QMC model was proposed to describe “new” physics beyond the standard 
hadronic picture in nuclear matter and hnite nuclei, we address the crucial question if the 
model is consistent with established results. As a hrst step we investigate the connection 
between the MQMC model and QHD. We demonstrate that in nuclear matter the MQMC 
energy functional is formally equivalent to the expression obtained in QHD with a general 
nonlinear scalar potential. The quark substructure is entirely contained in the scalar poten¬ 
tial and, in principle, this provides a tool to generate hadronic potentials based on a quark 
model. The explicit form of the hadronic potential depends on the model which is employed 
for the bag constant. Thus determining the parameters on the quark level is equivalent 
to calibrating the potential on the hadronic level which is a well established procedure in 
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nuclear matter calculations 16,17 


We apply our model to symmetric nuclear matter and compare the predictions with QHD 
mean-held calculations. We employ a version of QHD which contains quartic and cubic 
scalar self-interactions and which is calibrated to reproduce the same nuclear matter 
properties as the MQMC model. We show that different parametrizations and models for 
the bag constant lead to equivalent nuclear properties at low and moderate densities. In this 
region the MQMC model predictions are in excellent agreement with QHD. In contrast, the 
original version of the QMC model leads to substantially different results. Hence we are lead 
to the satisfactory conclusion that our generalized MQMC model predicts nuclear matter 
properties of the same quality as other established mean-held models. The key to this success 
is the correlation between the bag constant and nuclear matter properties. Conhrming the 
results of Refs. ||TI| , p!^ , we hnd that the bag constant has to be signihcantly smaller than 
its free-space value to reproduce the desired ehective nucleon mass at equilibrium. More 
importantly, the explicit value of the bag constant at equilibrium is model independent, i.e. 
independent of the details of the parametrization. 

The exact density dependence of the bag constant is essentially unknown. However, based 
on general theoretical arguments scaling relations among in-medium quantities have been 
proposed in the recent literature |]I^^. Since our approach establishes a direct connection 
between nuclear phenomenology and the quark substructure we investigate to what extent 
the proposed scaling relation between the bag constant and the effective nucleon mass is 
realized; at this point we encounter a clear model dependence. Although all our models 
are calibrated to produce identical nuclear matter properties we can build models which 
exactly follow the proposed scaling relation but also models which differ considerably. A 
similar model dependence can be observed for the bag radius. In accordance with previous 
studies [pT| , p!^ we hnd a sizable increase of the bag radius as a consequence of the decreasing 
bag constant. However, the quantitative predictions exhibit a strong model and parameter 
dependence. 

The outline of this paper is as follows: In Sec. H, we give a short description of the QMC 
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model and and snmmarize the relations which determine the nnclear matter properties. 
Section III is devoted to discnss the connection between the QMC model and QHD. We 
also discnss the calibration procednre. In Sec. IV, we apply onr model to symmetric nnclear 
matter. We compare onr resnlts with QHD and with the original version of the QMC model. 
Section V contains a short snmmary and onr conclusions. 


II. THE QUARK-MESON COUPLING MODEL 


In this section, we briefly summarize the relations which determine the nuclear equation 
of state in the quark-meson coupling model. For further details we refer the reader to 
Refs. @,0,0. 

In the QMC model the nucleon in the nuclear medium is described as a static, spherical 
MIT bag in which quarks couple to meson mean fields. In symmetric matter they are taken 
to be neutral scalar (a) and vector helds iV^). 

The energy of a bag consisting of three quarks in the ground state can be expressed as 

Ebag = 3 ^ - - + . ( 1 ) 

^ R R 3 ^ ’ 

where the parameter Z accounts for the zero point motion and B is the bag constant. The 
coupling of the quarks to the scalar held is inherent in the quantities Qq and x which are 
given by 


fig = -f {Rm*y 


,,, fn,-RmlY'\ , ^ 


( 2 ) 


and where m* = rrig — g^a denotes the effective quark mass and is the current quark 
mass. For simplicity we work in the chiral limit, i.e. mPg = 0. 

After correcting for the spurious center of mass motion, the effective mass of a nucleon 
bag is given by P] 


A2; = vMWsWp 


( 3 ) 
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For a fixed meson mean-field configuration the bag radius R is determined by the equilibrium 
condition for the nucleon bag in the medium 


dM*^ 

dR 


0 . 


( 4 ) 


In free space M^r can be hxed at its experimental value 939 MeV and the condition Eq. 
to determine the parameters B = Bq and Z = Zq. For our choice, Rq = 0.6 fm, the result 
for bI^"^ and Zq are 188.1 MeV and 2.03, respectively. 

The total energy density of nuclear matter can be written as 

^ dkP{k^ + + g^VoPN. (5) 

In symmetric matter (7 = 4) the Fermi momentum of the nucleons is related to the conserved 
baryon density by 


PN - —,k^ . 


( 6 ) 


The meson mean helds are determined by the general thermodynamic condition that they 
should make the energy per nucleon stationary. For the time-like component of the vector 
held this leads to the relation 


T/ 

g^rVo — —5-Pat , 

whereas the scalar held is determined by the self-consistency equation 


( 7 ) 


a = 


dk 


( 8 ) 


ml 27r^ JO {k? -|- 

The details of the quark substructure are entirely contained in the ehective coupling Cq{a 
which is related to the ehective nucleon mass by 

dM*^ 


cya) = 


da ’ 


(9) 


and which depends on the explicit form of the bag constant. Details can be found in 

Refs. 0 JT 3 . 
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III. GENERATING NONLINEAR MEAN-FIELD MODELS FROM MODIFIED 


QUARK-MESON COUPLING MODELS 


In the original version of the QMC model the bag parameters B and Z were held 

hxed at their free space values B = Bq, Z = Zq. The bag constant B is a nonuniversal 
quantity associated with the QCD trace anomaly. In the nuclear environment it is expected 
to decrease with increasing density as argued in Ref. [^]. At present, however, no reliable 
information on the medium dependence of B is available on the level of QCD calculations. 
On the other hand, much effort has recently been devoted to the study of effective models 
which approximate low energy QCD. The guiding principle for constructing such effective 
models are the symmetries and the symmetry breaking patterns of QCD, in particular chiral 
symmetry breaking. Typical representatives are Nambu and Jona-Lasinio (NJL) models 
and related chiral meson lagrangians [^. In this framework attempts have been made 


to model the QCD trace anomaly by introducing a scalar glueball held ||^. The concept 
of a bag constant arises naturally in these models and it is a common feature to predict 
a decreasing value of B when the density (or temperature) of the nuclear environment is 
increased p3[ . 


To account for this physics in the QMC approach two different models for the bag 
constant have been proposed [^,|^. A direct coupling between the bag constant and the 
scalar mean held 

,4 a 


B 


^- 9 . 


B 


( 10 ) 


S Mm- 

and a scaling model which relates the bag constant in the medium directly to the ehective 
nucleon mass 


B 




M, 


N 


( 11 ) 


Thus, rather than focus on the calculation of the bag constant in the medium the idea of 
the modihed QMC model is to parametrize the density dependence in terms of in-medium 
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quantities. A similar approach is used to construct nonlinear mean-field models in quan¬ 
tum hadrodynamics (QHD) where the unknown density dependence of the nuclear energy 
functional is parametrized by nonliner meson-meson interactions |]T6| , p^ . This task is cer¬ 
tainly more difficult on the quark level. The direct coupling model is inspired by NJL type 
nontopolgical soliton models for the nucleon |^, where a scalar soliton field is responsible 
for the binding of three quarks to form a nucleon. On the other hand, the scaling model is 
related to the idea of a general scaling relation between in-medium quantities as proposed 
by Brown and Rho lUJli. 

Generally, the parameter Z may also be modified in the medium |^. However, in 
contrast to the bag constant, there is less physical intuition how the medium dependence 
of Z can be cast in a model. Here we assume that the density dependence of Z can be 
disregarded and take Z = Zq. 

The original QMC model with B = Bq has a serious shortcoming. The predicted 
mean fields are much smaller than obtained in established relativistic mean-field models 
30^1- As a consequence, the effective nucleon mass is much too big |§J§]. Recently, 


it has been demonstrated that this shortcoming can be significantly corrected and that rel¬ 
ativistic nuclear-phenomenology can be recovered from the models given by Eq. m and 
Eq. (0) [0,|l^. However, the systematics of the predicted nuclear matter properties is still 
not satisfactory. The MQMC model produces acceptable values for the effective nucleon 
mass but values for the compressibility which are too high and vice versa 00. In our 
work we generalize the direct coupling model and the scaling model demonstrating that the 
resulting improved MQMC predicts nuclear matter properties of the same quality as in other 
successful relativistic mean-field models []. 

Furthermore, we will show that the relation between the QMC model and QHD is more 


^Experience has shown that an accurate reproduction of nuclear matter properties leads to realistic 
results when the calculations are extended to finite nuclei 116,27-30 
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general and direct than previously thought This will allow us a more consistent 

comparison between these different approaches. 

We start the analysis with a comparison of the energy density in Eq. (|^) with the corre¬ 
sponding expression in QHD [|^ 

Sqhd = ^ fj dk + u.(d,) , (12) 

where the effective nucleon mass is given by = Mat — gsfj)- The standard form of the 
nonlinear scalar potential is 





( 13 ) 


Although the cubic and quartic terms provide sufficient flexibility for an accurate calibra¬ 
tion of the nuclear equation of state []^, we will assume a general functional form for the 
potential. The scalar held is determined by the self-consistency equation 


“ 9i 


nk-p 


d(j) 


dk 




27 r 2 Jo {k^ + Mpy/^ 


(14) 


The equivalence of the two sets of equations, Eqs. (I) and (H), and Eqs. (|^) and (|TJ), can 
be demonstrated by performing a redehnition of the scalar held in the QMC model 


g„4,(a) = Mn - M-„{a) = M„ - ^. 


( 16 ) 


The transformation does not depend explicitly on the density if we assume 


as suggested in Eqs. ( 0 ) and ( 0 ). The coupling go is chosen to normalize the new held 
according to 


M + o{ay , 

cr^O 


and it is given by 


go = 




da 


(7=0 


( 16 ) 
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The transformation is well defined provided that the effective mass changes monotonically 
with the scalar held. Inverting the relation Eq. (|^ leads to a nonlinear scalar potential in 
the energy density Eq. 

1 


m,a (4) = U,(4) . 


(17) 


At this point the energy density in the QMC model is identical to the QHD expression 
Eq. (^) with a general scalar potential. For hnite nnclei the transformation also introdnces 
a change in the kinetic energy of the scalar meson. This will mainly lead to a different 
description of the snrface energy, bnt this effect is expected to be small. 

As shown in Ref. 0 a signihcant simplihcation arises for models with B = B{a) if 
= 0, i.e. if there is no direct conpling of the qnarks to the scalar held. In this case one 
obtains the scaling relation [|^ 


B 


M: 


N 


M, 


N 


(18) 


Thns, if the bag constant is of the specihc form 




9b 


a 




N 


(19) 


the transformation Eq. (|^ is linear leading to the original Walecka model with a 
qnadratic potential 

1 




( 20 ) 


The details of the qnark snbstructnre are entirely contained in the nonlinear potential 
Us{4>). This implies that if the bag constant, more generally the bag parameters, was known 
as a fnnction of the scalar held and the ehective mass, the steps leading to Eq. (0) wonld 
permit the prediction of potentials for hadronic mean-held models. 

On the other hand starting with a known potential, Eq. (0) dehnes the inverse held 


transformation to Eq. (|^ 


a(0) = 


217s(0) 




( 21 ) 
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After substituting this result in Eq. (|^, the QHD expression for the energy is exactly of 
the same form as in the MQMC model. The self-consistency condition Eq. (|^ changes to 


a = 


ml 27r^ Jo 
where the coupling is given by 


dk 


(A:2 + Mi^2)i/2 ’ 


( 22 ) 


a/2I/s 

C[a) = gsms-^ 

d4> 


(23) 


Thus, the MQMC model is formally equivalent to a nuclear mean held model with a field 
dependent scalar-nucleon coupling. Field dependent nucleon-meson couplings are often inter¬ 
preted as an indication for the compositeness of the nucleon [|6|,p!^ but from a modern point 
of view they appear naturally in hadronic models as a result of held redehnitions. In our case 
the held dependent coupling is equivalent to a subset of nonlinear meson self-interactions |^. 

Using the relation between MQMC and QHD one can determine the bag constant as 
a function of a. This approach is useful since it establishes a direct connection between 
nuclear phenomenology and the medium dependence of the bag parameters. 


To be more specihc let us assume B = B(a) and = 0. The scaling relation Eq. (h8|) 
together with the expression for the ehective mass in QHD 


then leads to 


= Mat - gs4>{cr) 


^ = (i_ 

Bq I Mat 


(24) 


(25) 


In general, the transformation which relates the two scalar helds is not known, but Eq. 
motivates the ansatz 


B 


= (l - with F(0) = 1 , 


(26) 


For a discussion of the equivalence between field dependent meson-nucleon couplings and non¬ 


linear meson self-interactions in the framework of QHD see Ref. []h| 
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which includes the original form Eq. ( 0 ). 

Similarly, the scaling model can be generalized by 

with G(l) = l. (27) 

We studied polynomial and Fade parametrizations for the unknown functions F and G 
which all lead to qualitatively similar results at low and moderate densities. In the following 
we will present results obtained by using a simple polynomial 

F(ct) = 1 + acr + /3(T^ , (28) 


for the generalization of the direct coupling model, and a Fade form 

g(m;) = ^ 


(29) 


l-a-6-c + aM*^ + ’ 

for the generalization of the scaling model. For notational convenience we will refer to the 


model in Eq. ([2^) as MQMCa and to the model in Eq. (|27|) as MQMCb. 

Since the parameters Bq and Z are hxed to reproduce the nucleon mass in the vacuum 
our models contain six free parameters. The parametrization of the bag constant contains 
the parameter k and the three couplings {gs, a, (3) and (a, b, c) for MQMCa and for MQMCb, 
respectively; in addition values for the quark-meson couplings and for the ratio g ^/are 
needed ^ Four of the six parameters can be chosen to reproduce the equilibrium properties 
of symmetric nuclear matter, which we take as the equilibrium density and binding energy 
(p^, —Co), the nucleon effective mass at equilibrium and the compression modulus 

{Kq). The hrst three of these are tightly constrained 0 , whereas the latter is not. The set of 
equilibrium properties used here H are listed in Table these are motivated by successful 
descriptions of bulk and single-particle nuclear properties |1T|J5(]|JHT[| . Since there are more 


^Strictly speaking, in nuclear matter only the ratios {gs!g%, (3j 9 %“^, Im^) are rel¬ 

evant for the calibration procedure. In order to vary the scalar coupling the mass of the scalar 
meson needs to be fixed. We choose = 550 MeV. 
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free parameters than constraints, we proceed as follows. We choose values for the coupling 

and for k and determine the remaining couplings by requiring that they reproduce the 
desired equilibrium properties. This is achieved by solving a set of transcendental equations 
that relate the parameters directly to the nuclear matter properties; this is a well established 
procedure in nuclear matter calculations | TB| , p!7 |. The original version of the QMC model 
with B = const, contains only two parameters, g'^ and g^/Tn^, which are usually chosen 
to reproduce the binding energy and the density at nuclear matter equilibrium. In our 
discussion we use the parameters given in Ref. |P]. For comparison we employ the QHD 
model with the nonliner potential in Eq.(|^. The parameters are determined to reproduce 
the same equilibrium properties as in the QMC model. 

Our primary goal is to study the influence of the in-medium bag constant on the equation 
of state (EOS) of nuclear matter. The freedom of varying the parameter k will allow us 
to investigate scaling relations between the bag constant and the effective nucleon mass. 
Moreover, models with g^ = 0 and 5^^ 7^ 0 describe different physical situations. For g'^ = 0 
the scalar meson couples only to the surface of the bag. Apart from an overall shift in the 
single-particle energies, due to the vector-meson mean held, the properties of the quarks 
inside the bag are not changed by the nuclear medium. In contrast, the quarks acquire a 
density dependent effective mass m* = — g^a for nonvanishing couplings g^. 

A decreasing effective nucleon mass in the original QMC model with B = Bq requires 
a positive quark-meson coupling {g^ > 0) which implies a negative effective quark mass. A 
similar effect can be observed in soliton models ||^. Close to the origin of the soliton the 
attractive scalar potential leads to a negative effective quark mass. But here this is less 
disturbing since the effective mass depends on the space coordinates and becomes positive 
in the outer region of the soliton. In the MQMC it is possible to chose g'^ > 0 or g^ < 0 and 
to generate positive or negative effective quark masses. However, we observe that the sign 
of the effective quark mass has no impact on the properties of nuclear matter. Models with 
either sign are qualitatively equivalent. 

A priori we have no specihc guidance on the allowed values of k and g^. We observe that 
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not all possible choices permit the reproduction of the desired equilibrium properties. For 
example, MQMCb has no solution for < 0.8. We analyzed MQMCa for 3 < k < 7 and 
0 < ~ 2, and MQMCb for 3 ^ ^ 5 and 0.8 ^ ^ 1.5. 

IV. NUCLEAR MATTER PROPERTIES 

The density dependence of the nucleon mass is controlled by the effective coupling Cq{a) 
in the self-consistency equation Eq. (H). For MQMCa this quantity is indicated in Fig. |l|. 
To normalize the curves at the origin we divided Cq{a) by g^ which was dehned in Eq. (|T6|). 
The effective coupling for QHD, given by Eq. (^31) , and for the original QMC model are also 
indicated. We emphasize that all parametrizations except the original QMC model reproduce 
the same equilibrium properties listed in Table |. For small values of the scalar held MQMCa 
and QHD are in good agreement predicting nearly constant effective couplings. In contrast, 
the curve for the QMC model steadily decreases and eventually becomes negative. 

The consequences for the effective nucleon mass can be studied in Fig. As expected 
from Fig. the masses for MQMCa and for QHD are nearly identical up to 1.5 nuclear 
matter densities. For B = Bq the effective mass decreases very slowly and, when the point 
Cq{a) = 0 is reached, increases again. The curves terminate at some maximum density, 
which corresponds to a maximum value of the scalar held. For g'^ = 0 this occurs at 
B = = 0. At that point the radius diverges. The curves for g^ ^ 0 terminate when x, 

the solution of Eq. (^, approaches zero. 

Fig.^ shows the ehective nucleon mass for MQMCb. Similar to what we observed for 
MQMCa, the curves are nearly indistinguishable from the QHD result at low and moderate 
densities. At higher densities MQMCb leads to a slowly decreasing ehective mass which 
approaches a nonzero asymptotic value. 

Fig.^ indicates the predicted nonliner potentials in Eq. (PT]) as a function of the trans¬ 
formed scalar held given by Eq. (|T5|). The value go(j)/Mj^ = 0.4 corresponds to the saturation 
point of nuclear matter. For gQ(l)/Mj^ < 0.6 the predicted MQMC potentials are almost iden- 
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tical to the nonlinear QHD potential in Eq. (|T^. In contrast the potential in the QMC model 


is much smaller. Here the saturation point is at go(j)/M]^ = 0.11. 

Due to the calibration procedure, the very good agreement of the different models and 
parametrizations at low and moderate densities is certainly not surprising. As clearly visible 
in Figs. the predictions of MQMCa and MQMCb vary considerably from QHD above 
the saturation point of nuclear matter. We also observe that different model types and 
different parametrizations within a specihc model are no longer equivalent at high densities. 
This is certainly an indication that an extrapolation into the regime of high densities might 
be problematic. 

In Fig. 1^ we show the binding energy curves of symmetric matter for MQMCa- The 
stiffness of the EOS is controlled by the parameters k and g^. The curves become softer for 
smaller values of k and stiffer if the quark-meson coupling is increased |1II|,p!^. Also indicated 
in Fig. ^ is the binding energy curve of the original QMC model. Although the compression 
modulus is only slightly lower [Kq = 223MeV) than in the other models {Kq = 250MeV) 
the EOS in the original QMC model is substantially softer. To understand this difference, 
it is useful to split the binding energy according to 

£ U U 

eo = — -Ma = — + ^ + ^, (30) 

Pn Pn Pn Pn 

with 


£^ = 


7 

2ti^ 


ckp 


'0 

and t/g = -mgCr^ . 


dkk^{k^ + - MtvPa , 


In Fig. ^ the three different contributions are indicated as a function of the density. Satura¬ 
tion arises as a competition between the decreasing effective mass which lowers the kinetic 
energy and the increase of the energy due to the increasing density. As expected, the curves 
for QHD and MQMCa are almost identical. The quantitative difference to the original QMC 
model is striking. To achieve saturation at a very high effective mass [M'^ = 0.89) only a 
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small contribution of the vector part U^, i.e. a small coupling is required. At high 

densities the term dominates and a smaller coupling leads to a softer EOS. This also 


explains the tendency of the QMC model to produce small vector mean helds . 

The in-medium bag constant as a function of the effective nucleon mass is indicated 
in Fig. 1^. We consider MQMCa in part (a) and MQMCb in part (b). To reproduce the 
desired effective nucleon mass at the saturation point, B is required to decrease substantially 
below its free-space value. The values of B/Bq are between the two lines and 

{M^/Mn)^. As mentioned in the last section, in the MQMCa model a special situation 


arises for = 0. In that case the bag constant scales as {M^/M^sy for all values of k |]T2 


We emphasize that all curves produce identical nuclear matter properties. Hence there is 
no compelling evidence from nuclear phenomenology that the bag constant should scale with 
a specihc power k, e.g., like proposed in Ref. |^. The tendency of the MQMC 

approach is to predict values of B which are slightly higher. This becomes more apparent 
in Fig. H where the bag constant at the saturation point as a function of g'^ is indicated. We 
examine MQMCa for three different values of the nucleon mass at equilibrium. The value 

= 0 correspond to the scaling behavior B/Bq = (M^/AAv)^. Because the calibration 
procedure determines the value of B, the curves are model independent, i.e. they depend 
only on the coupling but they are independent of the functional form of B = B{a). The 
curves in Fig. are also independent of the compressibility which merely determines the 
derivatives of B with respect to the scalar held in MQMCa and the derivative with respect 
to the effective nucleon mass in MQMCb. The key quantity here really is the effective 
nucleon mass. 

At this point, some caveats concerning the comparison with other approaches must 
be added. As mentioned earlier, the prediction of a reduced bag constant in the nuclear 
environment is also a common feature in effective models for low energy QCD, e.g. chiral 
models. Since none of these models can strictly be derived from QCD it is not clear to what 
extent results can be compared. Within the MQMC approach, it was one of our goals to 
demonstrate that it is possible to obtain nuclear matter results which are of the same quality 
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as in any other snccessful hadronic mean-held model. In contrast, chiral models are more 
concerned with the description of the underlying physics, i.e. the breaking and restoration 
of chiral symmetry. To the best of our knowledge, it has not been demonstrated whether it 
is possible to describe nuclear saturation and properties of hnite nuclei on that level. 

The corresponding bag radius is indicated in Fig. ^ As a consequence of the rapidly 
decreasing bag constant the MQMC models predict a picture of a signihcantly ’’swollen” 
nucleus ||ll|,[I^. This effect is most drastic in MQMCa for vanishing quark-meson couplings. 
As mentioned above, the effective nucleon mass and therefore also the bag constant van¬ 
ishes at some hnite value of the density. Consequently, the radius diverges at this point. 
Although all the models produce identical nuclear matter properties near equilibrium they 
yield signihcantly diherent radii, even at low densities. Also shown in Fig. ^ is the curve 
Rc = (3/47rp7v)^^^ which indicates the ’’critical” radius where the individual bags start over¬ 
lapping, signaling the breakdown of the simple bag model. In our model this occurs slightly 
above the saturation density. This behavior is in sharp contrast to the original QMC model 
which predicts a nearly constant radius. 

We emphasize that it is not possible to achieve reasonable nuclear matter properties and 
a small radius at the same time because the calibration procedure determines the value of 
the bag constant. As indicated in Fig. H, the acceptable range of M^/Mn = 0.6 —0.7 always 
requires a small value of B which in turn leads to a large radius. This is certainly a dilemma. 
On one hand we have demonstrated that the MQMC model can be improved and calibrated 
so that nuclear matter properties are accurately reproduced. The most important quantity 


here is the effective nucleon mass which is tightly constrained by nuclear observables [30 


On the other hand we found that the required bag constant in the medium leads to radii 
which are unreasonably high. Furthermore the radii are very sensitive to the model features 
and parametrizations. The size of a nucleon in matter is certainly a subtle concept and one 
might argue that physical observables do not depend on the bag radius [Q. However, it is 
an important phenomenological quantity in many nuclear physics issues where bag models 
are employed to describe features which depend on the intrinsic structure of the nucleon ( 
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see for example p4| ). 

Moreover, one motivation for generating the nuclear equation of state on basis of a quark 
model is the hope that such an approach is well suited to describe systems at high densities 
where quark degrees of freedom become important. However, because of the large radii 
and the picture of overlapping bags it is not clear how far the equation of state can be 
extrapolated into the high-density regime. 

V. SUMMARY 

In this paper we study properties of nuclear matter based on an improved quark-meson 
coupling model. This model describes nucleons as nonoverlapping MIT bags interacting 
through scalar and vector mean helds. Of central importance is the bag constant which we 
assume to depend on the density of the nuclear environment. We study two types of models 
for the bag constant in which the density dependence is parametrized in terms of in-medium 
quantities. In one we assume that the bag constant depends on the scalar held only and in 
another the density dependence is related to the effective nucleon mass. 

By performing a redehnition of the scalar held we demonstrate that the resulting energy 
functional corresponds to a QHD-type hadronic mean-held model with a general nonlinear 
scalar potential. In principle, this connection can be used to generate hadronic potentials 
from quark models. Moreover, a direct relation between nuclear phenomenology and the 
quark picture arises. We use this connection to motivate our models for the density depen¬ 
dence of the bag constant. 

For the explicit calculations we employ a polynomial and a Fade form to model the 
medium dependence of the bag constant. The unknown parameters can then be ht to 
properties of nuclear matter near equilibrium that are known to be characteristic of the 
observed bulk and single-particle properties of nuclei. 

Our basic goal is to study properties of nuclear matter. We investigate whether the 
models for the bag constant lead to results which are consistent with established hadronic 
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models. This is relevant in view of the hope to apply qnark models to describe “new” physics 
which goes beyond the hadronic pictnre. Becanse of the relation between MQMC and QHD 
we compare onr resnlts with a QHD model calibrated to prodnce the same eqnilibrinm 
properties. As the basic resnlt we hnd an excellent agreement between onr rehned MQMC 
models and QHD at low and moderate densities. In particnlar it is possible to reprodnce 
the desired saturation properties and large scalar and vector potentials as demanded by 
nuclear phenomenology. The central quantity here is the effective nucleon mass. Its accurate 
reproduction requires a rapidly decreasing bag constant below nuclear matter densities which 
leads to a satisfactory saturation mechanism. In contrast, the original QMC model predicts 
a very high effective nucleon mass and saturation is achieved on much smaller energy scales. 

Moreover, the calibration procedure determines the value of the bag constant as a func¬ 
tion of the effective nucleon mass in a model independent manner. However, this does not 
imply that the density dependence of B is well constrained. On the contrary, we hnd that 
models which reproduce identical properties of nuclear matter can be generated from dif¬ 
ferent parametrizations and models for the bag constant. In particular there is no clear 
evidence for a specihc scaling behavior of B. 

We observe a similar model dependence in the predicted bag radius. As a consequence of 
the decreasing bag constant the size of the nucleon signihcantly increases in all our models. 
This effect is even more drastic than in earlier calculations |Tl| , [I^ . We hnd that the nucleon 
bags start overlapping at densities slightly above the saturation point signaling the break 
down of the model. Moreover, diherent models and parametrizations which are equivalent 
near equilibrium produce diherent high-density equations of states. A similar picture arises 
in hadronic mean-held models a. Thus, even if quark degrees of freedom are incorporated 
it is unclear how far the nuclear equation of state can be extrapolated into the high-density 
regime. 

In recent applications the original QMC model was used to estimate the density depen¬ 
dence of quark condensates and other hadron masses [Q. In view of our results we expect 
that the predictions for these quantities change drastically if the density dependence of the 
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bag constant is taken into account. It is also important to investigate the uncertainties in 
these quantities which arise from the model and parameter dependence of the bag constant. 

In summary we conclude that an improved QMC model provides a very satisfactory 
description of nuclear matter near equilibrium. Phenomenological information in terms of 
equilibrium properties of nuclear matter can be used to constrain the unknown density 
dependence of the bag constant. On the hadronic level different models predict equiva¬ 
lent equation of states at low and moderate densities. This picture changes if microscopic 
quantities on the quark level are considered. Here a clear model and parameter depen¬ 
dence emerges. The difficulties rely on the detailed density dependence of the bag constant. 
Although these details have minor impact on the nuclear matter properties, microscopic 
quantities like the nucleon radius appear to be very sensitive. It is important to have these 
uncertainties under control before one can make reliable statements about the physics be¬ 
yond the standard hadronic picture. To achieve this goal more detailed empirical information 
on effective hadron masses and also additional observables from hnite nuclei might provide 
useful constrains to reduce model and parameter dependence in the future. 
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TABLES 


TABLE 1. Equilibrium Properties of Nuclear Matter 


{kp)° 

P% 


eo 

Ko 

1.3 fm“^ 

0.1484 fm-3 

0.60 

— 16.1 MeV 

250 MeV 
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FIGURE CAPTIONS 


FIG. 1. Normalized effective coupling Cq{a) for the model MQMCa, calculated with k = 4 
and for various values of g^. The corresponding quantity for QHD and for the original QMC model 
{B = const.) is also indicated. 

FIG. 2. Effective nucleon mass as a function of the density for the model MQMCa- Results 
for various values of k and are compared with QHD and the original QMC model {B = const.). 

FIG. 3. Effective nucleon mass as a function of the density for the model MQMCb- Results 
for different parameter sets are compared with QHD. 

EIG. 4. Predicted nonlinear scalar potential as a function of the transformed scalar field 
gQ(j) = Mtv — M^. The dotted and dashed curves correspond to MQMCa with g^ = 0 and g$ = I 
respectively. The dotted-dashed curves indicate the result for MQMCb- In addition the result for 
the original QMC model (B = const.) and the QHD potential is also indicated. 

EIG. 5. Binding energy as a function of the density for MQMCa- The dotted and dashed 
curves are calculated for k = 3,4, 5 from the bottom to the top. We also show the binding energy 
for the original QMC model {B = const.) and for QHD. 

EIG. 6. Kinetic energy S ^/vector potential U^/pn and scalar potential Us/P n a. function 
of the density. Results for MQMCa, QHD and the original QMC model {B = const.) are shown. 
The parameters for MQMCa are k = 3 and g'/ = 1. 

EIG. 7. Bag constant as a function of the effective nucleon mass calculated for k = A and for 
different quark-meson couplings g$. In part (a) we show the results for MQMCa- Here the curve 
{M/j/Mjq)^ corresponds to g% = 0. Part (B) indicates the results for MQMCb- 

PIG. 8. Bag constant at the saturation point for MQMCa as a function of the quark-meson 
coupling g%. The individual curves correspond to different values of the nucleon mass at equilibrium. 
The bag constant at equilibrium does not depend on k and the compressibility. 
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FIG. 9. Bag radius as a function of the density. At the critical radius Rc = (S/dvrpAr)^/^ the 


individual nucleon bags start overlapping. 
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